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ABSTRACT 

Standard cosmological models based on general relativity (GR) with dark energy 
predict that the Universe underwent a transition from decelerating to accelerating 
expansion at a moderate redshift z acc ~ 0.7. Clearly, it is of great interest to directly 
measure this transition in a mo del-independent way, without the assumption that GR 
is the correct theory of gravity. We explore to what extent supernova (SN) luminosity 
distance measurements provide evidence for such a transition: we show that, contrary 
to intuition, the well-known “turnover” in the SN distance residuals A ft relative to 
an empty (Milne) model does not give firm evidence for such a transition within the 
redshift range spanned by SN data. The observed turnover in that diagram is pre¬ 
dominantly due to the negative curvature in the Milne model, not the deceleration 
predicted by A cold dark matter and relatives. We show that there are several advan¬ 
tages in plotting distance residuals against a flat, non-accelerating model (w = —1/3), 
and also remapping the z—axis to u = ln(l + z); we outline a number of useful and 
intuitive properties of this presentation. We conclude that there are significant com¬ 
plementarities between SNe and baryon acoustic oscillations (BAOs): SNe offer high 
precision at low redshifts and give good constraints on the net amount of acceleration 
since z ~ 0.7, but are weak at constraining z acc ; while radial BAO measurements are 
probably superior for placing direct constraints on z acc . 

Key words: cosmological parameters - cosmology: observations - dark energy - 
distance scale. 


1 INTRODUCTION 

The A cold dark matter (ACDM) model has become well 
established as the standard model of cosmology, due to 
its very impressive fit to a variety o f cosmological obser¬ 
vations, including CMB a nisotropy llHinshaw et all 120131 : 
IPlanck Collaboration! 120141 '). large-scale galaxy clustering 
including the baryon acoustic oscillation (BAO) feature 
ll Anderson et all 12014!) , and the Hubb le diagram for distant 
supernovae (SNe; Betoule et ~aj|2014lb In ACDM and close 
relatives, the mass-energy content of the Universe under¬ 
went a transition from matter domination to dark energy 
domination in the recent past at a redshift z me ~ 0.33; the 
transition from decelerating to accelerating expansion, here¬ 
after z a cc, was somewhat earlier, at a redshif t z acc ~ 0.67. 
In ACDM, these are given by 1 + z acc = \/2Ul^/Ul Yn and 

1T Zme = ■y/flA/flm, so 1 + Zacc = \/2 (1 + Zme) . We see later 
that the value of z acc is relatively insensitive to dark energy 
properties, assuming standard GR and simple parametriza- 
tions of the dark energy equation of state. 

The most direct evidence for recent accelerated 


* E-mail: w.j.sutherland@qmul.ac.uk 


expansion comes from the many observations of dis¬ 
tant SNe at 0.02 < z < 1 .5; the early S N re¬ 
sults in 1998 dRiess et a jj 1 19981 : IPerlmutter et all il999f) 
began a rapid acceptance of dark energy, due also 
to p revious in direct ev ide nce f rom large -scale struc¬ 
ture (|Efstathiou. Suth erland fc Maddoxl Il990l ). the cluster 
baryon fraction (IWhite et aJ 1*19931 ') and the Hubble con¬ 
stant dFerrarese et all 1 19961 ). Strong independent support 
came from observation of the first CMB acoustic peak defin¬ 
ing a near-flat universe (Ide Bernardis et~a3l2000l : iBalbi et~all 
I2000l ). combined with decisive evidence for a low v alue o f 
f r om the 2dF Galaxy Redshift Survey (IPeacock et all 
l200ll : IPercival et all l2002h . In the past decade there has 
been a rapid improvement in the precision of observations 
in all these areas (see references above), most recently 
from the Planck, Baryon Oscillation Spectroscopic Survey 
(BOSS) and Supernova Legacy Survey (SNLS) projects. 
Current joint constraints are impressively consistent with 
ACDM with ~ 0.30 and Ho ~ 68.3kms -1 Mpc -1 
llAnderson et alll2014l : iBetoule et alll2014h . 


Many deductions in cosmology are based on six, seven 
or eight-parameter fits of extended ACDM to observational 
data, which generally show good consistency with the six- 
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parameter model and place upper limits on the additional 
parameters. However, given our substantial ignorance of the 
nature of dark energy, it is clearly interesting to ask what we 
can deduce with fewer assumptions, e.g. keeping the cosmo¬ 
logical principle while dropping the assumption of standard 
gravity. In particular, fitting models of GR with dark energy 
to the data produces a reasonably sharp prediction for the 
value of Za.cc', however, if the apparent cosmic a cceleration is 
due t o another cause suc h as modified g ravity (IClifton et all 
l2012h . a giant local void llCelerieij [20 071 or other, this may 
not necessarily hold; therefore, it is of considerable interest 
to see what constraints we can place on Zacc without assum¬ 
ing specific models._ 

It has been shown by e.g. IShapiro fe Turner] l|2006l l that 
the SN brightness/redshift relation does provide evidence 
for accelerated expansion independent of GR; but direct 
evidence for past deceleration is less secure. A number of 
other authors have explored GR-independent constraints 
on the cosmic exp ansion h istory , dark ener gy evolution 
and/or Zacc', e.g. ISahni fe Starobinskvl (12006 ) provide a 
broad review mainly foc used on dark energy reconstruction; 
ICattoen fe Visseil (l2008l l exp lore various distance de finitions 
related to z or y = 2/(1 + zk lCunha fe Limal (120081 ') derived 
constraints on a acc from SNe assuming simple parametriza- 
tions of deceleration parameter q{z)\ [Clarkson fe Zunckell 
(l2010h provide a method for non-parametric reconstruc- 
tion o f w(z) (mainly from future high-quality data); 
iMortsell fc Clarkson 1200911 prov ide non-parametric esti¬ 
mates of H(z)', and Nesseris fe Garcia-Bellidol d2013h give a 
comparison of several methods for estimating w(z ) from SNe 
data. Our work is partly related to these, but focusing more 
on the possibility of non-parametric constraints specifically 
on Zacc', where we overlap we are generally in agreement. 

The plan of the paper is as follows: in Section [2] we 
discuss the value of Zacc and the SN Hubble diagram, and 
the cause of the downturn in the latter. In Section[3]we point 
out several advantages of comparing SN residuals relative 
to a flat non-accelerating model. We discuss some future 
prospects in Section [I] and we summarize our conclusions 
in Section [5] Our default model is ACDM with S2 m = 0.300; 
Ho generally cancels except where stated. 


2 RELATION BETWEEN LUMINOSITY 
DISTANCES AND z acc 

2.1 The expected value of Zacc 

Here we note that the value Q of Zacc is now constrained 
rather well in flat wCDM models with constant dark energy 
equation of state w; for this model family, Zacc depends on 
only and w, and is given by 

1 + a a cc = [(-1 - 3iu)(l - B m )/n m ]“ 1/3u ' (1) 

1 In highly non-standard models, it is not guaranteed that z acc 
(defined by d = 0) is single-valued; e.g. if there were short-period 
low-amplitude oscillations in a, or a past accelerating phase tran¬ 
sitioned back to deceleration at a very low redshift, then in prin¬ 
ciple z aC c may be multi-valued. These possibilities appear im¬ 
probable and hard to test observationally, so we assume z acc is 
single-v alued (after th e CMB era) for the remainder of this paper; 
see also lLinderj (120 ] (I ) . 


(e.g. iTurner &~R icss 200^). This is shown in a contour plot 
in Fig.[l] It is interesting that in the neighbourhood of ~ 
0.3, w ~ — 1, the contours of constant z aC c are nearly vertical, 
thus Zacc is nearly independent of w and is well approximated 
by 

Zacc ~ 0.671 - 2.65 (Qm - 0.3) . (2) 

Qualitatively, this occurs because as w increases above — 1, 
there is less negative pressure hence less acceleration per unit 
Pde, but larger w gives higher pDE in the past; these effects 
happen to cancel (largely coincidentally) near the concor¬ 
dance model, so z acc is rather insensitive to w. This has 
positive and negative consequences: on the one hand, mea¬ 
suring Zacc is not useful for constraining w, on the other 
hand, the range 0.60 ^ z acc < 0.75 appears to be a robust 
prediction of wCDM, so if future data (e.g. direct measure¬ 
ments of H(z) from BAOs or cosmic chronometers, or new 
more precise SN data) were to empirically measure z acc out¬ 
side this range, it could essentially falsify the whole class 
of wCDM models. (Models with time-varying w such as the 
common model w{a ) = wo + w a (l — a) allow a wider range 
of Zacc, but these generally require z acc < 1 unless w a is dra¬ 
matically negative, w 0 < -1, which is disfavoured in most 
quintessence-type models). 

In Fig. |B we also show contours of (1 + Zacc)/E(z acc ), 
which is equivalent to the “net speedup” or integrated ac¬ 
celeration between z a cc and today; this is discussed later in 

m 

2.2 SN data 

For comparison with models, we use th e “Un ion 2.1” com - 
pilation of type-la SN distance moduli (ISuzuki et alll20f^ ~l. 
which contains 580 SNe of good quality spanning the range 
0.01 < z < 1.6. For plotting purposes we divide the sample 
into bins of approximately equal width in ln(l + z), while 
adjusting bin widths so that each bin contains ^ 20 SNe ex¬ 
cept at the highest redshifts; then, the mean distance mod¬ 
ulus residual and weighted average redshift are computed 
for each bin. The resulting binned data points are shown as 
‘Union 2.1’ in subsequent figures. 

We show a fit of this data set to flat wCDM models 
(with and constant w as the fit parameters; results 
of this fit are shown in Fig. [2] with a best-fitting point 
near = 0.28, w = —1.01. This shows the well-known 
degeneracy track between and w; here we note that 
the long axis of the track is quite similar to the contour 
(1 + Zacc)/E(zacc) ~ 1.15 in Fig. [T] this is discussed in later 
sections. 

We note th at a more recent SN la compilation has been 
produced bv lBetoule et all (|2014l ~) which includes more inter¬ 
mediate redshift SNe, more detailed photometric calibration 
and expanded treatment of systematic errors; however, the 
best-fitting parameters from the latter paper are within la 
of those above, so the slight difference is not important for 
the remainder of this paper. 

2.3 Fiducial models and A p 

The observations of Type la SNe are sensitive to the stan¬ 
dard luminosity distance Dl(z) for each SN, plus some scat¬ 
ter due to the intrinsic dispersion in absolute magnitude per 
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Figure 1. A contour plot of the acceleration redshift z acc , and (1 + z acc ) / E(z acc ), as functions of Q m ,w for flat wCDM models. The 
dotted horizontal line shows w = — 1. 

The solid black contours show z acc , in linear steps of 0.1 from 0.35 (right) to 0.95 (left). The dashed green contours show (l-\-z aC c )/E{ ■Zacc) 
(i.e. total net speed-up) in linear steps of 0.05 from 1.05 (upper right) to 1.35 (lower left). Selected contours are labelled. 


SN. In practice, the distant z > 0.1 SNe are compared to a 
local sample “in the Hubble flow” typically at z ~ 0.02 to 
0.05; for the local sample, peculiar velocities are assumed 
to be relatively small compared to the cosmological red- 
shift, so the value of Ho cancels with the (unknown) char¬ 
acteristic luminosity L c of a standardized SN. Thus, quasi¬ 
local SNe really constrain the degenerate combination h 2 L c 
or equivalently M c + 51og 10 h; and comparison of distant 
and local SN samples actually constrains the distance ratio 
Dl(z)/Dl(z ~ 0.03), rather than the absolute distance. 

The value of Dl{z) spans a very wide range over the 
redshift interval covered by SNe: from z ~ 0.03 to z ~ 2 
is a factor of ss 118 in distance or 10.3 magnitudes, while 
the differences between models are relatively modest: e.g. 
15 percent differences between ACDM and a zero—A open 
model, down to differences ~ 2 percent between ACDM and 
a w = —0.9 model. This implies that plotting Dl(z ) versus 
z directly is not very informative since model differences are 
very small compared to the plot range; therefore it is com¬ 
mon to present SN results as residuals relative to some fidu¬ 
cial model; residuals are often presented in distance modulus 
or magnitude units, i.e. 

A n(z) = 5log 10 (3) 

E>L,&d[Z) 

where is the value for some fiducial model. The choice 

of fiducial model is essentially arbitrary (up to small binning 
effects second-order in bin size); however, this choice can 
have a strong effect on the shape of the results and intuitive 
deductions, as shown below. 

One obvious choice of fiducial is ACDM itself; how¬ 
ever, this makes observed residuals (almost) flat-line, which 
does not translate readily into inferences on decelera¬ 
tion or acceleration. Another common choice of fiducial 
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Figure 2. The allowed region in the (ff m ,ru) plane from fitting 
flat constant-u) models to the Union 2.1 SN sample. Contours 
show the values of Ax 2 = 2.3, 6.0,10.6, corresponding to 68, 95 
and 99.8 percent confidence regions. 


model is the empty or Milne model, with D m = 0, 
Da = 0, fl k = 1, as used by many notable pap e rs e.g . 


Wa = U, iZ k = 1, as used by many notable pap e rs e.g . 
Riess et al (i 19981) ; i Leibundgufl d200il l: iR.iess et al ( 2004 ): 


Goobar fc Leibundeutl (|201ll ). The zero matter density 


means this is clearly not a viable model for the real Uni¬ 
verse, but it is a convenient fiducial model for two reasons: 


(i) It has a very simple analytic form for Dl(z), given by 

- iH 1 +1) ■ < 4 > 
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Figure 3. Distance modulus residuals relative to the Millie model for various cosmological models. The solid black lines show ACDM 
with f2 m = 0.27 (upper) and 0.30 (lower). Long-dashed red lines show the corresponding DO models (Equation 1121 with deceleration 
artificially turned off above 2acc- Dashed green lines show Friedmann models of historical interest: from top to bottom, a pure-vacuum 
model (Oa = 1); an open model with f2 m = 0.27, Da = 0; and an Einstein-de Sitter model (D m = 1). Dotted blue lines show constant-!/ 
models with q c = —0.6, —0.4, —0.33 respectively (top to bottom). 


hereafter we define A he to be distance modulus residuals 
relative to this. 

(ii) For a given Ho, the Milne model has the maximum 
luminosity distance among all Friedmann models with zero 
dark energy (assuming non-negative matter density). There¬ 
fore, observational evidence for distance ratios larger than 
the Milne model (positive A/rs) at any redshift is direct ev¬ 
idence that we do not live in a Friedmann model with zero 
dark energy. 

However, using the Milne model as fiducial has some draw¬ 
backs which we discuss in the next subsection; we sug¬ 
gest^ an improve d fiducial model in Section [3] (see also 
iMortsell fc Clarksonll2009l ). 


2.4 Downturn in distance residuals 

It is very well known that observed SN distance residuals 
are all significantly positive at 0.2 < 2 < 0.6, in agreement 
with the ACDM accelerating expansion. It is also fairly well 
known that ACDM models exhibit a turning point (a maxi¬ 
mum) in the A/te(«) relation. Fig.[3]shows that this turning 
point, hereafter Zt p , occurs at z ~ 0.50 for the = 0.300 
concordance model, and the predicted residuals then de¬ 
cline to a zero-crossing at z ~ 1.26. It is seen in Fig. [3] 
that the actual supernova data do hint at the existence of a 
turnover, with the three data points at z > 0.9 all slightly 
low compared to their predecessors. The actual evidence for 
this turnover is not decisive, but it is clearly somewhat pre¬ 
ferred by the data. The turnover occurs quite close to the 
theoretical transition epoch z a cc » 0.67, and it is therefore 
widely believed (at least anecdotally) that supernovae have 
directly detected the predicted cosmic deceleration at z > 1. 
We discuss some prior claims to this effect in Appendix m 

We demonstrate in the next subsection that the latter 
conclusion does not follow; specifically, while a downturn in 
A/ze is favoured by the data, the downturn predicted by 
ACDM is mostly caused by the negative space curvature 
in the fiducial Milne model, and cosmic deceleration makes 
only a minority contribution to the downturn. The fairly 
close match between ztp and z ac c is found to be largely co¬ 
incidental. 
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2.5 Cause of the turnover in A/ze 

Assuming homogeneity, the luminosity distance Dl(z) is 
given by 


(5) 

with E(z) = H(z)/Ho, and the function Sk{x) = 
sins;, x, sinhi for k = + 1 , 0 , —1 respectively, where k is the 
sign of the curvature (opposite to the sign of Dk, in the usual 
convention where fik = 1 — fltot)- 

It is convenient to factorize this so that 


D l {z) 

= (1 + z)Dr(z) 

(6) 


c r dz' r dz' 


Dr(z) 

H 0 Jo E 0') C Jo H(z') 

(7) 

X 

- HoDr{z) Vw 

(8) 


where Dr{z) is the comoving radial distance to redshift z\ 
and x is the dimensionless ratio between Dr(z) and the 
cosmic curva ture radius, which in a Friedmann model is 
R c = c/H q \/|flk | ■ We note that these distance results are 
still valid in a homogeneous and isotropic non-GR model, as 
long as the Robertson-Walker metric applies and we define 
f2k from the curvature radius via fik = A{HqR c / c)~ 2 , which 
is then not necessarily equal to 1 — f l to t- 

Looking at equation ©, the first (1 + z) factor is 
parameter-independent and due to time-dilation and loss of 
photon energy; these each give one power of (1 + z) -1 in flux, 
hence combine to (1 + z) in equivalent distance. The param¬ 
eter dependence of Dr(z) then factorizes into two parts, the 
D r (z) term dependent only on expansion history, and the 
factor Sk(x)/x which depends mainly on curvature and also 
(more weakly) on expansion history; this is asymptotically 
1 — kx 2 /6 for i < 1, or 1 + QkZ 2 /6 for z < 1. The factor¬ 
ization above is helpful to understand the relative impor¬ 
tance of curvature versus acceleration/deceleration on the 
distances and distance ratios. In the non-flat ACDM model, 
the combination of Planck+ BAO data requires |Dfc| < 0.008 
at 95 percent confid encfl (see equations 68a and b of 
iPlanck Collaboration! (I2014l 'l'l. which implies that the cur¬ 
vature factor is within 0.2 percent of 1 at the redshift range 
a < 1.5 of current SNe. 

It is now interesting to compare terms in equation © 
for the ACDM and empty models. In the case of the empty 
model, D r {z) = ( c/Hq ) ln(l + z), flk = +1, so equation © 
becomes 


Dl,e(z) 


(1 + z)—ln(l + z) 


sinh(ln(l + z)) 
ln(l + z) 


(9) 


2 We note that in non-GR models the standard limits on 
does not apply; however, if the true cosmology were a curved 
non-GR model, if \flk\ ^ 0-05 we would then require a rather 
close cancellation between curvature and non-GR effects in order 
to make the non-flat ACDM fits turn out so close to fi*. = 0. 
If we discard this possibility as an unnatural conspiracy, it is 
reasonable to assume \&>k\ < 0.05, and in that case the curvature 
factor Sk{%)/x « 1 d= 0.01 for jz < 1.5 for reasonable expansion 
histories. 


which easily simplifies to equation ©. However, it is more 
informative to keep the longer form of equation m since the 
rightmost fraction is a pure curvature effect; it is well ap¬ 
proximated by l + (ln(l + 2)) 2 /6 at z < 1. We next show that 
this term, not the transition to deceleration, is the dominant 
cause of the downturn in A/ze for models similar to ACDM. 

Considering the distance modulus residual A/z for any 
flat model relative to the empty model, we then have 


A he{z) 


5 logi 0 


Jo E(z') dz ' ln(l + z) 
ln(l + z) sinh(ln(l + z)) 


( 10 ) 


A uh[z) - A Hk{z) 


( 11 ) 


where we have broken the A /is into two additive terms, 
A hh(z) = 5log 10 [f* (1/ E(z'))dz' / ln(l + z )\ due to expan¬ 
sion histories, and Afik(z) = 5 log 10 [sinh(ln(l+a))/ln(l+z)] 
is the term due to curvature in the empty model (here de¬ 
fined so A/ifc is positive, thus it is subtracted in equation HU 
above). 

For illustration, we evaluate each of these terms for 
ACDM (with D m = 0.30) at two specific redshifts: we choose 
z a = 0.50 close to the turning point, and Zb = 1.26 to be 
the downward zero-crossing where A/ze(z) = 0. We then 
find A/ze( 0.50) = 0.1231 = 0.1822 - 0.0592 where the lat¬ 
ter two are Ahh and AfXk respectively. At zi = 1-26 we 
find —0.0005 = 0.2350 — 0.2355. Note that Ahh grows from 
z = 0.50 to z = 1.26, since although the expansion is de¬ 
celerating over most of this interval, the expansion rate a 
remains smaller than the present-day value; see below. 

For comparison purposes, it is useful to evaluate how 
much the predicted deceleration contributes to Ahh- for 
this we define another model set, hereafter DO, which ex¬ 
actly matches ACDM back to z a cc but with deceleration 
artificially switched off (q = 0) at z > z a c c : specifically, we 
define model DO by 

H(z) = HaCDm(z) if Z < Za.cc (12) 

= Hacdm(z acc )(!+*)/(! + ^acc ) if z > •2-acc 


The DO models are somewhat artificial, but have a continu¬ 
ous q(z) and are useful to isolate the relative contribution of 
ACDM deceleration on the observables. Also, they represent 
in a sense the closest possible match to ACDM among all 
possible non-decelerating models, so they are an interesting 
target to attempt to exclude observationally. The DO model 
(for D m = 0.30) is identical to the corresponding ACDM at 
z a , and at Zb we find AfiH = 0.2464 (and Afik = 0.2355 
again). Therefore, the actual brightening effect attributable 
to deceleration in ACDM is just the difference in A^h be¬ 
tween ACDM and DO, which is only —0.011 mag. This is 
smaller by a factor of 20 than the curvature effect; so, the 
bottom line of this subsection is that at z = 1.26, 95% of this 
downturn is due to curvature in the empty fiducial model 
(or 90% if we divide by the value A he = 0.1231 mag at its 
maximum). Either way, it is clear that the open curvature 
in the Milne model greatly dominates over deceleration as 
the source of the downturn in A he- 
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3 AN IMPROVED FIDUCIAL MODEL 


3.1 The flat non-accelerating model 


We have argued above that the presentation of distance 
residuals from the Milne or empty model is potentially con¬ 
fusing, since it leads to a generic curvature-induced down¬ 
turn in the residuals at 0 > 0.5 which occurs independent 
of whether the expansion really decelerated prior to that 
epoch. In this section we look at an improved fiducial model 
and demonstrate several advantages. 

In particular, the above discussion suggests a natural 
fiducial model is one with a constant expansion rate (decel¬ 
eration parameter q(z) = 0 , and H(z) = Ho(l + z) at all 
redshifts, as for the Milne model), but simply setting cur¬ 
vature to zero (equivalent to striking out the sinh in the 
equations above). This is equivalent to a Friedmann model 
with fl m = 0, f Ide = 1 and w = —1/3; hereafter model N 
for short . (This reference model has been employed previ¬ 
ously by ISeikel fe Schwar3 ll2008l l and iMortsell fe Clarksonl 
(M), but appears to be rather uncommon in the liter¬ 
ature.) Again, this model is not realistic due to the zero 
matter density, but it is useful since it has both zero de¬ 
celeration and zero curvature. This model straightforwardly 
gives 

D l ,n{z) = -^-(1 + z) ln(l + z) . (13) 

n o 

We now define the distance ratio for any other model, yu(z), 
as the ratio Dl(z)/Dl,n(z), therefore 


Vd(z) 


HqDr(z) S k (x) 
cln(l + z) x 


(14) 


For an almost-flat model at a < 1.7 we can again neglect 
the curvature term as very close to 1 (as per footnote in 
Sect. 12.51 

Thus, for flat models the distance ratio becomes 


2/d (z) 


1 r dz' 

ln(l + z) J Q E@j 


(15) 


For many purposes below, it is more convenient to change 
the redshift variable to u = ln(l + z), which gives 


Vd(u) 


i ri + z' 

U Jo E ( z ') 


du ; 


(16) 


as usual u',z' are dummy integration variables, not deriva¬ 
tives, and 2/d(m) means 2/d (z = e“ — 1 ). 

This ?/d is directly related to A p,H above via A (j.h(z) = 
51og 10 yu(z), but several results below are simplified if we 
choose not to apply this log. Since i/d (a) is fairly close to 1 
in reasonable models, this is anyway rather close to a linear 
stretch A/t h ~ 2.17(i/d — 1). 

Since E(z)/( 1 + z) is just the expansion rate at a rela¬ 
tive to the present day, i.e. a{z)/a{z = 0 ), the integrand of 
equation m is just the inverse of this; i.e. yn(z) measures 
the average value of (a)o /a with respect to ln(l + «), over the 
interval from the source to the present. It is more convenient 
to work with averages of (1 + z)/E(z) rather than 1 /E(z), 
since the former varies much more slowly with redshift: for 
our default ACDM model, (1 +z)/E(z) reaches a maximum 
value of 1.153 at 2 acc — 0.67, crosses 1 again at 2 ~ 2.08, 
and declines to 0.895 at 2 = 3. 

Note also that since (1 + z)/E(z) contains the inverse 


of a, while 2 increases backwards in time, derivatives of 
(1 + z)/E(z) have the same sign as a, i.e. positive for 
acceleration. In fact the standard deceleration parameter 
q = —a/(aH 2 (a)) is given by 


g(u) = _i. ln fL+4 N | 

’ du \E(z)) 
which is useful below. 


(17) 


3.2 Useful properties of 2 /d 

The above definition of 2/d is simple and intuitive, and 
we show below that it enables a number of useful non- 
parametric deductions, as follows: 

(i) It is clear above that a value of yo(z) > 1 at any 
2 implies the past-average of a was less than the present 
value, i.e. acceleration has dominated over deceleration over 
this interval (note, this is not strictly the same as requiring 
a > 0 at the present day); this feature is similar to the Milne 
fiducial model above. 

(ii) It is easy to see that if q(z) is always negative over 
some interval 0 SC 2 ^ 2 i, then (1 + z)/E(z ) is a strictly 
increasing function of 2 , and therefore so is yn{z); i.e. a flat 
model which is non-decelerating at 0 < 2 < 21 cannot have a 
turnover in j/d at 2 ^ 2 i, regardless of the specific expansion 
history. The converse of this is that if a turnover in yn{z) is 
observed, this implies a transition to deceleration must have 
occurred within the interval, i.e. we can definitely conclude 
2 acc < Ztu independent of the functional form of E(z). Also, 
if a turnover exists at ztu, differentiating equation ca im¬ 
plies that the value of 1/a at ztu was equal to its average 
value (w.r.t. u ) across the interval from ztu to today. 

(iii) We can improve on the results above using the Mean 
Value Theorem: specifically, if had a known value yo(zi) = 
yi, this theorem implies that there exists some 2 < 2 i with 
(1+z)/E(z) > 2/1 ; i.e. the cosmic expansion rate has speeded 
up by at least a factor of j/i since some z < z i, independent 
of the functional form of E(z). For a more realistic case 
where we measure an average value of i/d in a finite bin, 
e.g. (//d) = y averaged between 21 < 2 < 22 , we can use 
the Mean Value Theorem twice, first, there exists some 2 m 
within this bin with yn{z m ) = y, and secondly there exists 
some 23 < 2 m < 22 satisfying (1+ 23 )/E( 23 ) > y. The above 
argument applies for exact knowledge of y, neglecting error 
bars; however, it is clear that the same argument also applies 
if we insert an observational lower bound for y. 

(iv) Also, it is interesting to ask a reverse question: if the 
expansion was decelerating at all 2 > 2 ac c, does this imply 
that a turnover in 2 /d(z) must exist ? The answer appears 
to be ‘almost always’: it is possible to build a contrived ex¬ 
pansion history where q(z ) crosses from negative to a small 
positive value, then asymptotes back to zero from above at 
high 2, so (1 +z)/E(z) tends to a constant from above; in this 
contrived case we can have deceleration at all 2 > 2 ac c while 
2 /d (z) monotonically increases to the same constant. How¬ 
ever, if we assume non-infinitesimal deceleration, q(z) ^ +t 
for all 2 > 21 and some positive value e, it is readily proved 
that 2 /d (a) must have a turnover at some 2 (though not nec¬ 
essarily in a readily observable range). 

(v) Differentiating equation m and rearranging gives 
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Figure 4. As Fig. [5] but with the horizontal axis now linear in u = ln(l + z). 



Figure 5. The distance ratio yu( z ) defined in equation cu for various cosmological models. As in Fig. [3] solid black lines are ACDM 
models with = 0.27 (upper) and 0.30 (lower). Long-dashed red lines are corresponding DO models, with deceleration artificially 
switched off. The short-dashed green lines are four Friedmann models of historical interest: from top to bottom, vacuum-dominated 
(Q m = 0 ,^a = 1); empty (Milne); open (Q m = 0.27, Dde — 0); and Einstein-de Sitter (£l m = 1)- Dotted blue lines are three constant-g 
models with q = —0.6, —0.4, —0.33 (top to bottom). Points with errorbars show the binned Union 2.1 SNe data. 
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l + z 
E(z) 


yv{u) + u 


dyo 

du 


(18) 


This gives us a direct graphical implication: taking the tan¬ 
gent to the curve of yn(u) at any point mi and extrapolating 
the tangent line to u = 2ui gives us directly the value of 
(1 + z)/E(z) at zi = exp(ui) — 1. 

Differentiating again shows that the transition to acceler¬ 
ation occurs when d 2 yu/du 2 = — (2/u)dyn/du; however, as 
is well known the need to take a second derivative of noisy 
data implies that this is not a very useful method for directly 
estimating u ac c- 

(vi) Substituting from equation C3 above leads to the 
compact results 


yo(u) 

q{u) 


u Jo 

2 d vp> 

du 


2/d( u ) + 


u d2y P 

d u. ~ 

dy-p 


(19) 

( 20 ) 


this shows that qo = —2(dyu/du)(0), but also that as u in¬ 
creases we get increasing weight from the second-derivative 
term, so it becomes increasingly more challenging to con¬ 
strain q(u) directly from numerical derivatives of data with 
realistic noise. Even for optimistic 1% error bars on j/d 
in bins Ait = 0.1, we get order-unity errors on d 2 yo/du 2 , 
so free-form reconstruction of q(u) is essentially impossible 
given realistic errors; the best we can do is assume some 
smooth few-parameter model for q(u) and fit. 

(vii) From equation m it clearly follows that for two 
measurements at redshifts corresponding to iti, U 2 we have 


U2yo(u2) — uiyo(ui) 


112 — 111 


112 — Hi 


l + z' 
E(z') 


du' 


( 21 ) 


where the right-hand side (RHS) is the average of (1 + 
z)/E(z) between the endpoints; therefore we can estimate 
this average as a linear combination of the two values at 
the ends; this is simple with respect to combination of error 
bars, and does not assume 112 — 111 is small. 

(viii) We now show another useful property of yn- for 
any flat model with 5 ( 2 ) = constant (of either sign), the 
second derivative d 2 yu/du 2 with respect to u is everywhere 
non-negative. For such a model, denoting q c as the constant 
value of q , we have H(z) = Ho(l + z) 1+9c . This easily leads 
to 


D l (z) 
2/d (2) 
2 /d (u) 


£( 1 + ,)A [ (1 + 

-1 (1 + z)~ 9c - 1 
q c ln(l + z) 

1 - e~ qcU 
q c u 


Now differentiating twice with respect to u gives 

d 2 yo = -1 e~ qcU {u 2 g 2 + 2 ug e + 2 ) — 2 ~ 

du 2 q c u 3 

2 |~ 2 — e~ p (p 2 + 2 p + 2 )~ 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


where we define p = q c u. The function in square brackets 
above is positive for all p, thus the above second derivative 
is everywhere non-negative for any value of q c with either 


sign, and is zero only if q c = 0 and yo = 1. For the cases 
of interest here, we are mainly interested in — 0.6 < q c < 0 
at 0 < 11 < 1 , hence — 0.6 < p < 0 ; the square-bracket term 
evaluates to 1/3 for p = 0 and 0.53 for p = —0.6, so for any 
reasonable q c model the second derivative is then between 
0.33 Qc an d 0.53g/, i.e. small, positive and slowly varying 
with u. 

This has a useful consequence: if q(u) were in fact any con¬ 
stant, then the graph of yrt(u) versus u must always show 
positive curvature (concave from above). Conversely, if the 
observed data points for yrt(u) exhibit significant negative 
curvature over some interval, we can conclude that q(u) in¬ 
creased with u at some point within the observed interval, 
again regardless of the specific functional form. (Note this 
does not necessarily imply that q(u) became positive, merely 
that it increased with u i.e. was less negative in the past.) 

We note that in the above points, items (i)-(iv) apply 
whether we choose z or u as the redshift variable, but items 
(v)-(viii) only apply with u as the variable; this suggests the 
latter is preferred. 

For an illustration of the current data, we plot yo(u) 
against u = ln(l + z) in Fig. [4] Although this is a sim¬ 
ple transformation of the x— axis from Fig. [5] the qualita¬ 
tive appearance is somewhat different due to the non-linear 
transformation, i.e. higher redshifts become squashed. The 
apparent “knee” in the ACDM models around 2 ~ 0.5 in 
Fig. 0 is significantly smoothed out with the u— axis, and 
both ACDM models now look very close to simple parabolas 
(see below). Also, the constant -5 models change curvature 
from negative in Fig. [5] to small and positive in Fig. U as 
derived above. Comparing to the data, it is clear that the 
SNe data points do marginally prefer a negative curvature 
in yn(u), but not overwhelmingly so. 

To quantify this, we fit three models to the yn(u) 
data points: a linear model, a quadratic, and the family of 
constant -5 models above; we find that the quadratic model 
is preferred over the linear model by Ay; 2 = 3.5 for 1 extra 
degree of freedom (d.o.f.), while the quadratic is preferred 
over the best constant -5 model by Ay 2 = 5.7 for 1 extra 
d.o.f. This indicates that negative curvature in t/p (increas¬ 
ing 5 ) is preferred, but only at around the 2 a significance 
level. We expand on the quadratic model below. 


3.3 A quadratic fitting function for (1 + z)/E(z) 

Here we note that it is interesting to consider a fitting func¬ 
tion where 1 /a is a quadratic function of u, specifically 

4^4 = 1 + biu - b 2 u 2 (27) 

E{z) 

with arbitrary constants 61 , 62 , and u = ln(l + 2 ) as before. 
The minus sign above is chosen so that positive 61 , 62 leads 
to recent acceleration and past deceleration as anticipated, 
with u a cc = 61/262 from equation (1171) . This fitting function 
is not physically motivated, but is useful since it provides 
a very good approximation to models similar to ACDM at 
u < 1, (2 < 1.72) (see Appendix [Clfor an approximate expla¬ 
nation of this property), and it gives several simple analytic 
results below. 

Fitting this function to the default ACDM (1 + z)/E(z) 
over 0 < u < 1 (2 < 1.72) gives best-fitting values 61 = 
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0.569, 62 = 0.530 with an rms error of 0.28 percent, and 
a worst-case error of —0.8 percent. (This fit becomes sig¬ 
nificantly worse above z > 2 , and has a catastrophic zero¬ 
crossing at u ~ 2 (z ~ 6.4), but it is good over the range 
accessible to medium-term SN data.) The functional form 
(EZ1) gives simple relations between M acc and the turnover in 
2/d; it easily gives 


2/d(u) 

= 1 + ifeiu - i b 2 u 2 ; 

( 28 ) 

q{u) 

— b i + 2fe 2 u 

1 + biu — 62M 2 ’ 

( 29 ) 

D l (u ) 

= -^(1 +z)(u+ ibiu 2 - 

^ 2 u 3 ) ; (30) 

so 2/d(w) 

is also an exact quadratic 

in this case. The q[u ) 


behaviour is approximately linear at moderate u, so this 
model is fairly similar to the model q(a) = ?o + q a (l — a) 
used elsewhere. Equation ESI) with values 6 i, 62 as above 
matches the exact numerical Dl{z) for ACDM with very 
high accuracy, a maximum error only 0.13 percent back to 
u — 1; this error is substantially smaller than for E(z), due 
to the integral for Dl- 

We find that the results above also work well for wCDM 
models in the region 0.2 < < 0.4,—1.2 < w < —0.8; 

thus, it is interesting (and partly a coincidence) that any 
wCDM model within the presently-favoured range leads to a 
2 /d (m) curve virtually indistinguishable from a quadratic, to 
around the 0.2 percent level i.e. comparable to the line thick¬ 
ness in Fig. [4] This gives another helpful feature: any proof 
of ‘percent-level’ deviation of yrt(u) from a simple quadratic 
would signify a failure of wCDM. 

We now look at the relation between 2 acc and the turn¬ 
ing point in 2 /d- In the above model equation ED with 
&i ,&2 > 0 , recall the acceleration epoch is u acc = 61/262 , 
hence (l+2 aC c)/-F(2 aC c) = l + 6f/462 ; while the maximum in 
2 /d occurs at Utp = 36i/ 4&2, at height yn(ut P ) = l+36f/1662. 
So, in this model 2 acc is directly related to the location ut P 
of the maximum, and (1 + 2 aC c)/E( 2 a cc) is directly related 
to its height, via 


^acc 

— 5 ^acc — (1 T ^tp) ^ 1 j 

(31) 

1 + -Zacc 
E ( Zacc ) 

= 1 + ^( 2 /D(utp) - 1 ) 

(32) 

without 

requiring to solve for 61 , 62 - 



This suggests that for other reasonably smooth 
parametrizations of E{z) such as wCDM models, we may 
expect equations m and ED to hold approximately, rather 
than exactly as above. In our default ACDM model, the ex¬ 
act values are 2 acc = 0.671, (1 + 2 a cc)/E(z acc ) = 1-1530, 
while from numerical evaluation of ut p and yo(ut P ) the 
RHS of the above equations evaluate to 0.693 and 1.1525 
respectively; thus equation EH is quite good, while equa¬ 
tion (13211 is an excellent approximation. More generally, we 
have tested these for wCDM models (constant w) with the 
results shown in Fig. [6] ; this shows that equation (13211 re¬ 
mains very accurate for a substantial range around the con¬ 
cordance model. 

We have also tested linear-;/ models q(a) = qo+q a {l—a), 
and find that equation ED is accurate to better than 0.01 for 
reasonable values of qo,q a , while equation m is somewhat 
worse but generally good to a few percent. For varying -w 



Figure 6. This figure shows the peak value of ?/d against 
the integrated acceleration (1 + z aC c)/E(za, C c), for a grid 
of wCDM models. The differing point styles show w = 
— 1.2,—1.1,—1.0,—0.9,—0.8 as indicated in the key. For each 
value of w we show seven points with Q m = 0.24, 0.26,..., 0.36 
in linear steps of 0.02; in each case these run from f2 m = 0.24 at 
upper-right to 0.36 at lower-left, so the central point is Q m = 0.30. 
The dotted line (not a fit) is equation (1321 1. 


models of the form w(a) = wo + w a (l — a), these approxi¬ 
mations remain good for w a ^ 0 but become somewhat less 
accurate for negative w a , especially for w a < —0.5. 

The summary here is that equation El is generally an 
excellent approximation for constant-ui models, and a good 
approximation for varying-w if w a is not too negative; while 
equation ED is fairly good at the few-percent level. 

These approximations are useful since the right-hand- 
side of equations ED and El are in principle directly ob¬ 
servable: it is clear from Fig. [4] that the location of the 
possible maximum in 2 /d is relatively poorly constrained, 
but if the suggestion of negative curvature in 2 /d is real and 
persists as expected to higher redshifts, then the SNe dat- 
apoints imply that yn(u) is probably approaching a maxi¬ 
mum value ~ 1.10 — 1.14 at ut P < 1; if so, this would give a 
direct and reasonably model-independent inference of the in¬ 
tegrated acceleration (1 + 2 a cc)/-E(z acc ) « 1.13— 1.18. This 
provides a useful intuitive explanation of the ridge-line of 
fl m versus w observed in Fig. [2] 

To summarize this subsection, we find that u>CDM 
models with constant w near the concordance model are 
very well approximated by the above fitting functions, i.e. 
very close to simple quadratics in 2/d(w), and thus equa¬ 
tions ED and ED provide quite accurate approximations 
relating the observable turning point in 2 /d to 2 acc and the 
net acceleration. 

Finally, in Appendix [B] we use the fitting function of 
equation ED to provide a simple and accurate ‘computer- 
free’ approximation to the luminosity distance in wCDM 
models. 
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% 

Figure 7. The allowed region in the (qo,q a ) plane from fitting 
models with q(a) = qo + q a (X — cl) to the Union 2.1 supernova 
data. Elliptical contours show the values of A% 2 = 2.3,6.0,10.6 
corresponding to 68, 95 and 99.8 percent confidence regions. The 
sloping lines bound the region of no deceleration and the region 
z acc < 2, with the wedge between these giving z acc > 2. The line 
along the major axis of the ellipse is illustrative and gives a pivot 
value q(a = 0.815) = —0.36 at a = 0.815 (z — 0.227). 


3.4 Linear q(a) models 

Here we briefly consider the two-parameter model family 
with deceleration parameter q given by a linear function of 
scale factor a, i.e. 

q{a) = q 0 + q a (l - a) (33) 

for constants qo, q a . This model has been used before by var - 
ious authors fe.g. ICunha fc Limall2008l . ISantos et allEoilf) . 
since it is simple, fairly flexible and can produce a fairly good 
approximation to the behaviour of many dark energy models 
at z < 2. We have fitted this parameter pair to the Union 2.1 
SN data, with best-fit values at (-0.62,-1-1.40) and the re¬ 
sulting likelihood contours shown in Fig.Q as expected, neg¬ 
ative qo is required at very high signih c ance. (This agrees 
well with a similar figure in ISantos et all (|201lf 0. The figure 
also shows lines bounding the regions of no past deceler¬ 
ation qo + q a < 0, and the region z acc < 2 equivalent to 
qo + 2q a /3 > 0; the wedge between these lines corresponds 
to a transition redshift z acc > 2. This plot shows that the 
no-deceleration region is disfavoured at around the 1.3 ct con¬ 
fidence level, but there is a region inside the wedge z aC c > 2 
which is allowed at around 0.8cr. In this wedge, no deceler¬ 
ation occurs within the redshift range of observed SNe, so 
the inference of deceleration relies on a linear extrapolation 
of the q(a) model beyond the range of SNe. This generally 
agrees with our previous conclusions, that a trend of less 
negative q at higher redshift is clearly preferred, but there 
is negligible evidence from SN data alone for an actual tran¬ 
sition to deceleration within the observed range. 


4 DISCUSSION 

It is instructive to blink back and forth between Figs [3] 
El El above: although from a parameter-fitting perspective 
there is no difference since the residuals (data-model) are 
all the same, from the perspective of visual intuition about 
expansion rate there are rather striking differences between 
these three Figures. Clearly, Fig. [3]shows a fairly convincing 
turnover in the data points; while in Fig. El the data shows 
negligible evidence for a turnover, but a reasonably con¬ 
vincing change in slope to a broad near-flat “plateau” above 
z > 0.6. Finally, in Fig. Elthe ACDM models are extremely 
close to parabolic (i.e. near-constant negative second deriva¬ 
tive), while the data points show near-linear behaviour with 
a reasonable but non-decisive indication of negative curva¬ 
ture; the constant-q models show weak positive curvature 
as derived earlier in equation As we argued earlier, 

the turnover in Fig. [3] is largely attributable to the negative 
space curvature in the Milne model, not due to actual decel¬ 
eration. Figs El and El show a much more gradual turnover in 
the ACDM models, while the DO models show the expected 
gradual rise; clearly the current data are completely unable 
to discriminate between ACDM and DO models. We suggest 
that Fig.|4]is the most informative due to the various useful 
intuitive properties outlined in § 13.21 above. 

The above conclusions seem somewhat unexpected: 
there is a widespread view (see Appendix lAl) that the SN 
data has convincingly verified the expected deceleration of 
the universe at z > 1. However from the discussion above, 
the SNe data are almost entirely inconclusive on the sign of 
q at z > 0.7, and even a constant-g model with q{z) ~ —0.4 
back to z > 1 is only excluded at the ~ 2.5cr level which is 
significant but not overwhelming. Thus, there is moderately 
good evidence for q increasing in the past, but concluding 
that q actually crossed zero to a positive value relies strongly 
on a smooth extrapolation of this trend, and is therefore 
model-dependent. 

Conversely, if we assume GR, almost all the acceptable 
models imply significant deceleration at z > 1. Essentially, if 
we assume GR with the weak energy condition and a value 
of >0.2, then the eightfold increase in p m back to z = 
1 combined with the much slower increase in dark energy 
guarantees matter domination and deceleration at z > 1; 
in this case deceleration at z > 1 is mainly a prediction 
of GR, rather than a feature directly required by data. For 
the value of z acc it is important to keep clear the distinction 
between an extrapolation based on GR parameter-fitting, or 
an actual detection purely based on data. 

It is clear that the CMB does provide much stronger 
constraints due to the long distance lever-arm: if we assume 
the standard sound horizon length inferred from Planck, 
then we deduce yo(z ~ 1090) ~ 0.44, which clearly requires 
a turnover and hence deceleration. However, since the CMB 
only gives us one integrated distance to z ~ 1090 spanning 
seven e—folds of expansion, while the supernova data con¬ 
strains only the last one e—fold of expansion, it would be 
straightforward to construct ‘designer’ expansion histories 
with some extra deceleration hidden in the un-observed six 
e—folds to offset an absence of deceleration back to z ~ 1.7. 
This is clearly contrived, but would not directly conflict with 
any available Dl (z) data. Therefore, even adopting the stan¬ 
dard distance constraint from the CMB, we do not yet have 
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a GR-independent proof that the expansion was actually de¬ 
celerating at 1 < z < 2; this is clearly the most probable and 
least contrived interpretation, but loopholes remain. 

We note that recent BAO results do provide significant 
evidence for decele ration; from t he firs t detection of BAOs in 
the Ly-a forest bv lBusca et all (l2013h and comp arison with 
lower-redshift measurements. [Busca et~afl (120131 1 quote 


E(z = 2.3)/3.3 
E(z = 0.5)/1.5 


1.17 ±0.05 


(34) 


which is a 3.4 a detection of deceleration between the above 
two redshifts (though this does assume an external WMAP7 
curvature constraint, which introduces some slight level of 
GR-dependence). However, the desirable goal of verifying 
that Zacc < 1 as expected is considerably more challenging, 
since the expected change in a between z = 0.67 and 1 is 
only 1.7 percen t in our default model. The Euclid spacecraft 
liLaureiis et alllioTl l is predicted to get sub-percent mea¬ 
surements of r s H(z) at a range of redshifts 0.9 < z < 1.8, 
which looks very promising for a direct model-independent 
result, while improved ground-based measurements span¬ 
ning 0.3 < z < 0.9 would also be highly desirable. 


5 CONCLUSIONS 

We summarize our conclusions as follows: 

(i) The predicted value of z acc is rather well constrained by 
current data within wCDM models, and is mainly sensitive 
to rather than w; this implies that a direct measurement 
of Zacc is not helpful for measuring w, but is potentially an 
interesting test of icCDM versus alternate models such as 
modified gravity. 

(ii) Contrary to intuition, the (probable) downturn in SN 
residuals relative to the empty Milne model does not provide 
convincing evidence for deceleration. The predicted down¬ 
turn is strongly dominated by the negative space curvature 
in the Milne model, and the actual deceleration in ACDM 
makes only a small minority contribution to the downturn. 

(iii) There are many advantages to presenting SNe dis¬ 
tance residuals relative to a flat coasting model (fl m = 0, 
fl DE = 1, w — —1/3), and also in changing the horizontal 
axis from z to « = ln(l ± z) as in Fig. [4] This presenta¬ 
tion enables a number of robust non-parametric deductions 
about expansion history based on the global shape of the ob¬ 
served residuals i/d (u), without needing specific numerical 
derivatives of data or fitting functions. Notably, a turnover 
in this plot is decisive evidence for deceleration, while any 
negative curvature in the data points is evidence for higher 
q in the past. 

(iv) If a turning point in yn(u) is observed, then we can 
infer z acc from its location and (1 + z a cc)/-E(z aC c) from its 
height from Eqs. (1311321) ; the latter relation holds to very 
good accuracy in the case of uiCDM models, slightly de¬ 
grading in the case of large negative w a ■ 

(v) For the case of znCDM models near the concordance 
range, the model curves of i/d (w) are remarkably close to 
simple quadratics to an rms accuracy < 0.3 percent, sig¬ 
nificantly better than present data. This provides a simple 
intuitive visual test for potential deviations from wCDM. 

(vi) For constraining expansion history, there are signif¬ 
icant complementarities between SNe and BAO (or cos¬ 


mic chronometers): the SNe have a precise local anchor 
at z ^ 0.05 and therefore place strong constraints on the 
integrated acceleration, e.g. giving robust lower bounds on 
the value of 1.7/ E(0.7) ^ 1.1. However, the combination of 
the integral in SNe distances and the broad maximum in 
(1 + z)/E(z) around the acceleration transition implies that 
SNe are weak at giving model-independent constraints on 
z a cc- In contrast, BAOs offer direct access to H(z) without 
differentiation and are therefore potentially stronger at con¬ 
straining z acc ; but they have limited precision due to cosmic 
variance at z < 0.25, and they are therefore weaker at con¬ 
straining the total integrated acceleration, most of which 
occurs at 0 < z < 0.3. 

It is clearly important to get a good cross-anchor between 
SN measurements and BAO measurements for c onstrai ning 
the absolute distance scale; as argued by e.g. ISutherlandl 
J2012ll . precision measurements of both SNe and BAO at 
matche d redshifts would be very useful for this; see also 
iBlake et all (l201lh for a slightly different but related ap¬ 
proach. 
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APPENDIX A: PREVIOUS CLAIMS OF 
DECELERATION 

Here, we provide a short discussion of previous claims con¬ 
cerning evidence of past deceleration from SN data; these 
are mostly in press releases or the semi-popular literature, 
but have had a significant influence. One of the earliest 
such claims appears to be a quote from A. Riess on NBC 
Newi0 , 2001 April 02, ‘the new supernova, dubbed SN 
1997ff, confirms that the universe began speeding up rel¬ 
atively recently’. A notable Scientific American article (Feb. 
2004) by Riess & Turner includes the quote ‘the observa¬ 
tions (six SNe > 7 Gyr old) confirmed the existence of an 
ear ly slow down period.’ This appears to be partly based on 
iTurner fc Riessl (j2002j ). and that paper fits two classes of 
model: first ruCDM models (in which case deceleration oc¬ 
curs almost ‘by assumption’ for reasonable values of fl m ); 
and secondly a two-parameter model in which q{z) follows a 
step transition between two constant values, an early value 
52 to a late-time value gi; also the transition redshift was 
artificially fixed at 2 = 0.4 or 0.6, so this model set is quite 
restrictive and not very represent ative of plausible dark en¬ 
ergy evolution. Also, the title of iRiess etali (120041 ) (R04) 
contains the phrase ‘Evidence for past deceleration...’; that 
paper is (as of 2014) the most-cited astrophysics paper pub¬ 
lished in 2004, and has thus been highly influential. Specif¬ 
ically, R04 Section 4.1 considers a two-parameter model 
q(z) = qo + z(dq/dz ) with constant dq/dz, and find that 
dq/dz is positive (implying past deceleration) at above the 
95% confidence level. Converting to z a cc = —go/ (dq/dz), 
R04 derived z acc = 0.46 ±0.13. However, we note that there 
are several possible caveats in this result: firstly a constant 
dq/dz model is somewhat unphysical since it leads to di¬ 
vergent q at large 2 ; more realistic models like ACDM have 
dq/dz decreasing with 2 , so a linear q(z) model tends to 
underestimate 2 acc ; a linea r q(a) relatio n as in Section 13.41 
ICunha fc Limal (l2008l l and ISantos et all ll201 if) is probably 
more realistic. Secondly, the choice of uniform priors in 
go, dq/dz leads to a prior density which is steeply rising to¬ 
wards small 2acc• If the true model is close to ACDM, both of 
these effec ts may te nd to p ull the 2 a cc estimate low. Thirdly, 
as seen in ICunha fc Limal (l2008l 'l, inclusion of more recent 
SNLS SN data also shifts the likelihood contours slightly 
towards smaller dq/dz\ most of their samples exhibit some 
non-decelerating regions inside the 95% confidence contour. 

From our discussions above, it appears that there has 
been a tendency to overstate the strength of evidence for 
actual past de c elerat i on in SN data; the r esults of R04, 
ICunha fc Limal (|2008l ), iMortsell fc Clarkson! (|2009h and this 
paper agree that there is reasonable evidence (~ 2 <j) that 
q(z) was less negative in the past than today, but the GR- 
independent evidence for an actual zero-crossing (i.e. tran¬ 
sition to deceleration) is relatively weak, and sensitive to 
the choice of parametric form for q{z). Thus improved data 
is highly desirable to prove past deceleration at high confi¬ 
dence. 


3 http://nbcnews.com/id/3077854 
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APPENDIX B: A SIMPLE APPROXIMATION 
FOR D l 

Here, we note that we can also use the fitting function equa¬ 
tion OH) to obtain a simple ‘computer-free’ approximation 
to Dl(z) for constant-ui wCDM models near the concor¬ 
dance model, which is remarkably accurate up to 2 < 1.7. 
The procedure goes as follows: 

(i) For given the standard Friedmann equation 

gives the value of 2 acc as in equation |l]), hence (1 + 
Zacc )/E(z a . cc ) follows. 

(ii) Given those two values above, we can then readily 
solve for the pair ( 61 , 62 ) in equation (1271) which reproduce 
the same position and value of the turning point in (1 + 
z)/E(z); the result is 61 = 2[(1 + 2 ac c)/E( 2 aC c) — l]/M aC c, 
62 = 6 i/ 2 u aC c where u acc =ln(l + s aC c). 

(iii) Finally inserting the above constants 61,62 in equa¬ 
tion gives our simplified approximation for Dl- 

Since this procedure matches only the turning point in 
(1 + z)/E(z), as expected it results in slightly different val¬ 
ues of 61 , 62 compared to the previous case in § 13.31 where 
we numerically fitted 61,62 to the Friedmann (1 + z)/E(z) 
function over the full range 0 < u < 1. Thus we get a less 
accurate approximation to Dl, but the accuracy still turns 
out surprisingly good for this back-of-envelope level approx¬ 
imation. For the case of the concordance model = 0.30, 
w = —1, the recipe above gives 61 = 0.5960, 62 = 0.5804, 
hence the approximation becomes 

D l {z) ~ {c/H 0 )(l + z)(u + 0.298u 2 - 0.1935 u 3 ) ; (Bl) 

comparing this to the quasi-exact numerical Dl gives an 
rms error of 0.14 percent and a worst-case error 0.32 per¬ 
cent across the range 0 < u < 1 (z < 1.72). The ac¬ 
curacy improves for w > — 1 and degrades for w < — 1 , 
but remains < 0.3 percent rms across the preferred ranges 
0.27 < fl m < 0.33, —1.2 < w < —0.8. The recipe above 
is substantially more accurate than a traditional third-order 
Taylor expansion in z, which rapidly becomes poor at 2 > 1. 


which corresponds to 2 < 0.86 and u < 0.62 in our default 
model. Also, differentiating equation dm) gives 


dq d 2 ( 1 + 2 \ 

dii ~ \E(z) ) ’ 


(C3) 


thus a slowly-varying dq/du leads to near-quadratic depen¬ 
dence for (1 + z)/E(z) versus it. We note that this is partly 
coincidental for parameters near the concordance model, 
since the present-day value qo ~ —0.55 is near the end of 
the timespan when +0.1 > q > —0.6. 


This paper has been typeset from a TgK/ ETj^X hie prepared 
by the author. 


APPENDIX C: RELATION BETWEEN q AND 
JERK 


Here, we provide a short argument why the simple fitting 
function of equation (1271) works surprisingly well at u < 1. 
With the deceleration parameter q defined as above, and the 
dimensionless jerk parameter j defined by 


. _ d 3 a/dt 3 
J = am {a) ’ 

it is shown by e.g. iBolotin et all :2012i i that 


(Cl) 


dq 

du 


j - q{2q + 1) 



(C2) 


(which is model-independent, assuming only that the deriva¬ 
tives exist). I 11 ACDM models, j = +1 independent of time 
(for negligible radiation content); this implies that dq/du 
had a maximum when q = — 1/4, and was slowly varying be¬ 
tween 7/8 and 9/8 over the period with —0.60 < q < +0.10, 
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